Mesoscopic "Rydberg" atom in a microwave field 
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We establish analogy between a microwave ionization of Rydberg atoms and a charge transport 
through a chaotic quantum dot induced by a monochromatic field in a regime with a potential barrier 
between dot contacts. We show that the quantum coherence leads to dynamical localization of 
electron excitation in number of photons absorbed inside the dot. The theory developed determines 
the dependence of localization length on dot and microwave parameters showing that the microwave 
power can switch the dot between metallic and insulating regimes. 
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The dynamical localization of quantum chaos is a 
generic physical phenomenon induced by quantum coher- 
ence which leads to suppression of diffusive wave packet 
spreading originated from dynamical chaos in the classi- 
cal limit (see e.g. and Refs. therein). It has been first 
seen in numerical simulations of the kicked rotator model 
and, later, of a more realistic system of excited hydro- 
gen atom in a microwave field The first observation 
of this phenomenon was achieved in microwave ionization 
experiments of hydrogen and Rydberg atoms 0, 0, IE 
while more recent experimental progress allowed to real- 
ize the original kicked rotator model with cold atoms in 
laser fields and detect with them the dynamical localiza- 
tion of chaos 

For a hydrogen atom in a microwave field the diffu- 
sive excitation in energy appears only above certain field 
threshold where the integrability of unperturbed motion 
is destroyed and the classical dynamics becomes chaotic 
0. However, in certain systems the internal dynamics 
can be fully chaotic on an energy surface so that clas- 
sically the diffusive excitation starts for arbitrary small 
fields. As an example of such systems we may quote 
complex molecules Rydberg atoms in a magnetic 

field [IH or chaotic Sinai billiards ^ n sucn systems 
the quantum eigenstates are chaotic and their level spac- 
ing statistics is described by the random matrix theory 

ES IH H 

The excitation of such quantum systems 
by a monochromatic field represents a diffusive process of 
one-photon absorption/emission transitions with a rate T 
which is given by the Fermi golden rule. The quantum 
interference effects lead to the dynamical localization of 
this diffusion with the localization length given by ^3 : 
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Here 1$ is measured in the number of photons of 
monochromatic field with frequency to, p c is the level 
density of unperturbed system and it is assumed that 
hiop c > 1 and 1$ > 1. The extensive numerical sim- 
ulations of a microwave ionization of internally chaotic 



Rydberg ato ms ( chaos is produced by static magnetic or 
electric field) [la] showed that the relation (JJJ works even 
in extreme regimes when up to a thousand of photons is 
needed to ionize one atom. 

The result is rather general and can be used not 
only for atoms but also for chaotic billiards in a mi- 
crowave field. In fact, during the last decade the con- 
ductance properties of such billiards realized with 2D 
electron gas quantum dots of micron size have been stud- 
ied in great detail both experimentally and theoretically 
(see reviews 0, and Refs. therein). In such ex- 
periments typical parameters correspond to a dot size 
a ~ 1pm and electron density n e ~ 4 • 10 11 cm -2 . With 
the spin degeneracy the level number at the Fermi energy 
is rip — n e A/2 sw 2000 where the dot area is A ~ a 2 . 
These values of a and n p correspond to those of a Ryd- 
berg atom with the principal number n ~ 100. Indeed, 
the atom size is a ~ n 2 a,B ~ 0.5/xm and the level num- 
ber, inside a set with fixed magnetic quantum number 
m z ~ 1, is riF ~ n 2 /2 ~ 5000 (see e.g.[l^j). Here as is 
the Bohr radius and m z is preserved for a linearly polar- 
ized field. Thus, a quantum dot with the above parame- 
ters can be viewed as a mesoscopic "Rydberg" atom. A 
tunneling barrier of large hight U between the dot and 
the outcome lead (see Fig. 1) creates an effective ion- 
ization potential similar to the one of Rydberg atom. A 
small barrier Ui n between the dot and income lead can be 
used to control electron confinement inside the dot. The 
microwave field will switch the dot conductance from in- 
sulating (Tiidlfj, <C U) to metallic (huil^ > U) regime. 

To understand in a better way the requirements for 
microwave parameters we note that for the 2D electron 
gas the Fermi momentum pp and energy Ep are given 
by Pp = 2Trn e h 2 and Ep = p F /2m with m = 0.067m e 
for nGaAs. The average level spacing in the dot is A = 
l/pc = Ep/np = 2ttH 2 /mA. For the typical values of a 
and n e used above we obtain Ep ~ 100K, A ~ 0.05i^ 
with the Fermi velocity vp — pp /m ~ 2-10 7 cto/s. Hence, 
the collision frequency is v c /2ir — vp/2a ~ 100GHz 
that is much larger than the frequency corresponding 



FIG. 1: (a) Dot billiard at A = 3/8 with the probability 
density of the eigenstate at the Fermi level uf = 2001 (see 
text); density is proportional to grayness. Dashed stripes 
show schematically income (left) and outcome (right) leads, 
(b) Sketch of tunneling barriers Ui n and U between income 
lead, dot and outcome lead . 

to one level spacing A/(2tt7j) ~ 1GHz. Thus for a 
microwave field with oj/27t < 100GHz ac-driving is in 
the classical adiabatic regime and according to the adi- 
abatic theorem the energy excitation is exponentially 
small for dot billiards with intcgrablc dynamics. How- 
ever, for dots with chaotic dynamics a microwave field 
with hw > A leads to diffusive excitation in energy E 
with the rate Dp = (SE) 2 /5t « Ti 2 uj 2 T. fn analogy with 
a chaotic atom in a linearly x-polarized field we have 
dE/dt — eujx cos tot and De ~ (eaco) 2 jv c where e is the 
field strength multiplied by the electron charge. This 
gives 

Z « 2ir X e 2 a 2 /(hv c A) « 16 X e 2 (-4/.4o) 5/2 KoK) 1/2 

(2) 

where x is a numerical constant and in the right part 
Aq = lum 2 , n eQ — 4 ■ 10 11 cm -2 and e is measured in 
V/cm. It is also convenient to write @ in the form 1$ ss 
x( ea / Ef) 2 ti 3 J 2 which shows that 1$ may be large even 
when ea <C Ep. 

To check the validity of the above estimates and ob- 
tain all numerical coefficients we study numerically the 
dynamical photonic localization in the billiard which is 
given by a conformal quadratic map of a unit circle pro- 
posed by Robnik [l^. Its shape in polar coordinates is 
r(</>) = R(l + 2Acos(/>) with a parameter A S [0,1/2]. 
For A > 1/4 the billiard is non-convex and almost fully 
chaotic, while for A = 1/2 it is rigorously known to be 
ergodic and belongs to the class of K-systems . We re- 
strict our studies mainly to the case with A = 3/8 shown 
in Fig. 1. For the quantum evolution we consider only 
even states of the billiard. These states are symmetric 
with respect to reflection transformation y to — y and 
a microwave field linearly polarized along x-axis gives 
transitions only inside this symmetry class. We take the 
Fermi level to be np — 2001 in this symmetry class (to- 
tal quantum number around 4000) that corresponds to 
Ep = Vp/A = 12586.2 where for numerical simulations 
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FIG. 2: (color online) Probability distribution (w n ) over the 
billiard eigenstates n after a pulse of 2000 microwave field 
periods (averaged over last 500 periods, n is counted from the 
Fermi level np). Here e = 1200, ui — 20 (red/black curve) and 
e = 300, ui = 160 (green/gray curve). Initially all probability 
is in the state n = shown in Fig. 1. The straight dotted lines 
show the theoretical exponential decay with the localization 
length ©. 

we use the usual billiard units with h = R = 2m = 1. 
Then near Ep the average level spacing for this symmetry 
class is A = 6.266. The numerical method introduced in 
[hH allows to find efficiently the billiard eigenstates and 
the dipole matrix elements between them (the eigenstate 
at the Fermi level is shown in Fig. 1). Therefore, the 
quantum evolution induced by a microwave field can be 
numerically integrated directly in the eigenbasis. This 
approach allows to follow the quantum excitation over 
thousands of microwave periods. The integration time 
step was set to At = 10~ 4 and we ensured that its mod- 
ification did not affect the excitation probabilities. In 
the numerical simulations the transitions below np were 
suppressed, since in a dot all states below np are occu- 
pied by electrons, and up to n tot — 2000 levels above 
were taken into account. The escape in outcome lead is 
modeled as absorption of all probability above the level 
number np + ni after a microwave pulse of finite du- 
ration. This corresponds to the tunneling barrier hight 
U = n/A between the dot and outcome lead. 

Examples of probability distribution w n over billiard 
eigenstates n after a long microwave pulse are shown 
in Fig. 2. They clearly show exponential localization 
of probability which can be approximately described by 
w n ~ exp(—2n/l)/l. The localization length in the basis 
n is related to the photonic length 1$ = lA/u). For large 
to = 160 the distribution shows a chain of equidistant 
peaks corresponding to one-photon transitions with the 
distance between peaks 8n = u>/A. For u> ~ A the peaks 
disappear and w n decays in a homogeneous way. 

The localization length 1$ can be expressed through 
the correlation function of dynamical motion inside the 
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FIG. 3: (color online) Dependence of the scaled localization 
length £</,/e 2 (left axis) on the microwave frequency oj. Numer- 
ical data for e = 1200, 1200, 1200, 600, 300 (from left to right) 
are shown by points with error bars, the full curve shows the 
theory J3J with the classical spectral density S(uj/uj c ) taken 
from |20j| for the billiard of Fig. 1 (sketch), here uj c » 224, 
A = 6.266. Right axis shows the scale for S(uj/uj c )- The 
dashed curve shows the classical spectral density S(uj/uj c ) for 
the 3-disks billiard (sketch) ; u) is adjusted to the same collision 
frequency uj c as for the full curve case. 



billiard. Indeed, as discussed in [20j in the semiclassical 
regime the dipole matrix elements of x can be expressed 
via the spectral density of dynamical variable x(t). Us- 
ing the definition of De and the expression for dE/dt we 
obtain De — e 2 uj 2 R 2 S (uj /uj c ) /2uj c . Here uj c = vf/R and 

S(k) — limT-^oo S-t dT^(r)e lKT /2T is the dimension- 
less spectral density of £(r) = x{t)/R with r = ui c t and 
k = u>/u> c . Together with JIJ this result leads to 
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S(uj/uj c ) 



(3) 



In the billiard units uj c = 2\fEp and for our value of Ep 
we have uj c w 224 and l^/e 2 ~ 5/448. The classical spec- 
tral density S(uj/uj c ) has been found numerically in |20| 
and we show it in Fig. 3. As it is seen, S(oj/oj c ) has one 
maximum with 5(1) « 1.8 and it goes to a constant value 
5(0) « 0.004 near zero. Using © we can determine the 
quantum localization length 1$ from the classical value of 
5. Without any fit parameters this relation gives a good 
agreement with the numerical data of Fig. 2. 

To check the theory Q in more detail we determine the 
value of I,/, from the fit of exponential probability decay: 
(w n ) oc exp(— 2nA/uol c f > ) (w n is integrated in one-photon 
interval to extract 1$). The dependence of l^/e 2 on u> is 
shown in Fig. 3 demonstrating a good agreement with 
the theoretical formula 10 |2j]- The dependence on e at 
a fixed uj is shown in Fig. 4. It is also in a good agreement 
with the theory while 3> 1 which is assumed by the 
semiclassical approach ;22j. 

The nontrivial property of Eq. © is that is directly 
dependent on the dynamical spectrum S(u>/u> c ). For the 
mesoscopic billiards discussed above uj <C uj c and it is 
interesting to have billiards with large value 5(0). Since 
5(0) is simply the integrated auto-correlation function 



1000 



100 - 



10 




0.01 



0.1 

£R/E F 



FIG. 4: (color online) Localization length I,/, as deter- 
mined from numerical data (points) versus dimensionless field 
strength €R/Ef for uj — 20. The straight line gives the theo- 
retical dependence 



of £(r), a possible choice is a billiard formed by three 
touching disks which exhibits a slow polynomial decay of 
correlations (see e.g. H^). In Fig. 3 we show the classical 
power spectrum S(u>/u) c ) versus u) for the 3-disks billiard 
(see sketch, the opening between touching circles has a 
length 0.05i?). In this case the spectrum has a plateau 
at small uj that gives the value of 5(0) by a factor 250 
larger compared to the billiard of Fig. 1. The comparison 
of with © gives x = 5(0)/4 (a = 2R, v c = ituj c /2). 
For the 3-disks billiard we have x ~ 0.25 and 1$ w 1 for 
e = 0.5V /cm and A — Aq, n e = n e Q. 

In fact the billiard shape is of primary importance for 
the efficiency of microwave excitation. Indeed, a decrease 
of A from A = 3/8 to A = 0.1 makes the unperturbed dy- 
namics inside the billiard quasi-integrable. In the latter 
case a microwave field acts in an adiabatic way and the 
current through the dot is reduced by orders of magni- 
tude compared to the chaotic billiard as it is shown in 
Fig. 5 24]. This figure also shows that for a chaotic 
billiard the current through the dot is strongly reduced 
as soon as HujI^ becomes smaller than the barrier hight 
U = rtjA. Therefore, a conductance g between leads 
can be efficiently changed by varying a microwave field 
power. The conductance appears due to diffusion in en- 
ergy which is rather similar to spatial diffusion in meso- 
scopic quasi-one-dimensional wires. Using analogy with 
the latter case where g = E c /A [23, we can write the 
microwave conductance of the billiard as 



g ~ Z exp(-2iV / /^)/iV / 



(4) 



where E c = HDe/U 2 is the Thouless energy for diffusion 
in energy space, Ni = U /huj is the number of photons 
required to pass over the lead barrier U and an effective 
level spacing between quasienergy levels is A e / ~ A/Nj. 
As in 25] it drops with the sample size Nj and g ~ 
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FIG. 5: Dependence of the current / through a dot billiard 
on a scaled microwave field strength eR/EF at u — 20, uf ~ 
2001. The current I, shown in arbitrary units, is determined 
as the probability on billiard eigenstates, within 2500 < n < 
4000, absorbed after a field pulse of 2000 microwave periods 
duration. Two curves correspond to fully chaotic billiard at 
A = 3/8 (top) and to quasi-integrable billiard at A = 0.1 
(bottom). The arrow marks the position where the absorption 
border is reached by localization m = 500 = huil^/A. 



E c /A e f. For lef, ^S> Nj the dot is in the metallic regime 
while for 1$ <C Ni it is insulating. Thus a microwave 
field can efficiently switch on and off the conductance of 
a chaotic quantum dot blocked by a tunneling barrier U 
between the dot and outcome lead. To distinguish clearly 
the dependence Q from the Arrhenius activation law at 
finite temperature T we need to have T <C hul^. 

Of course, above we used a rather simplified model ne- 
glecting interactions between excited electrons and finite 
temperature of a dot. These effects destroy quantum co- 
herence and dynamical localization and should be taken 
into account when comparing the theory with the experi- 
ment. However, we expect that a transition from metallic 
to insulating behavior induced by a microwave field (see 
(0J) is a robust phenomenon and thus it can be observed 
experimentally as it has happen with a microwave ioniza- 
tion of Rydberg atoms. Present experimental techniques 
allow to observe effects of microwave radiation on elec- 
tron transport in mesoscopic dots (see e.g. [IE 113 ) that 
makes possible experimental investigations of the effects 
discussed here. 

We thank Kvon Ze Don for discussions which origi- 
nated this work. 
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